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SUMMARY 


The reconstructed magnitude of radiance fluctuations (or fields) from dis- 
crete radiometric measurements is subject to errors caused by the aliased noise 
generated when the radiometer under samples the spatial detail of these fluctua- 
tions. Quantitative results indicate that aliased noise can be a large source 
of error in the reconstruction process when the sampling intervals equal the 
instantaneous field of view (IFOV) of the radiometer (i.e., for contiguous 
coverage) . This noise can be reduced about 2 orders of magnitude by shaping 
the effective IFOV (formed by the photosensor aperture and signal electronics) 
of a scanning radiometer and modestly decreasing the sampling intervals rela- 
tive to the IFOV (to provide some overlap between successive measurements) . 

In contrast, spatially smoothed (i.e. f low-pass filtered) reconstructions of 
the discrete measurements provide only a relatively modest decrease in this 
noise at the cost of spatial resolution. 


INTRODUCTION 

The reconstruction of spatial radiance fluctuations (or fields) is subject 
to errors caused by the aliased noise that is generated if spatial details have 
been undersampled. The resultant degradation of image quality obtained with 
line-scan devices (e.g., television cameras and optical-mechanical scanners) 
has received considerable attention in the literature (refs. 1 to 9) . These 
investigations have emphasized spatial resolution: the reproduction, detec- 

tion, recognition, and identification of detail. The accurate reconstruction 
of the radiance magnitude, or of spatial averages of the radiance magnitude, 
has been of little concern. The degradation of radiometric accuracy due to 
aliased noise has (to the authors 1 knowledge) never been quantitatively deter- 
mined, except for a few illustrative examples (refs. 6 and 10). 

This investigation formulates and computes the magnitude of the aliased 
noise that degrades the accuracy of continuous reconstructions of discrete 
radiometric measurements. It provides evaluations, in particular, of the mag- 
nitude of aliased noise as a function of the spatial (or angular) response and 
sampling interval of the radiometer and of the resolution of the reconstructed 
measurements. To produce generally useful results, a Wiener spectrum that is 
approximately representative of a wide range of scenes is used to characterize 
the (random) radiance fluctuations. This spectrum can be derived if the radi- 
ance fluctuations are assumed to be a random set of two-dimensional pulses 
whose width and magnitude obey the Poisson and the Gaussian probability density 
function, respectively. 


SYMBOLS 

a,b photosensor aperture shape functions, m (see fig. 3) 



vQ > 


g(x,y) 

(U,W) 

Ji 

L(x,y) 

L(U,W) 

m(x,y) 

m(x,y) 

m(x,y;X,Y) 

m(U,U);X,Y) 

m,n 

M,N 


spatial function 

spatial frequency spectrum of g(x,y), m 2 
first-order Bessel function 
radiance fluctuations, Win -2 

frequency spectrum of radiance fluctuations, W 
continuous signal, W 

continuous signal of smoothed measurements, W 
discrete measurements, W (see fig. 1) 
frequency spectrum of measurements, Wm 2 (see fig. 5) 
measurement count along x- and y-axis, respectively 
constants that determine rectangular dimensions of resel 


r = I (x-) - x 2 ) 2 + (y] - y2) 2 J 1//2 

R resolution cell (resel), m 2 

x,y measurement coordinates, m (see fig. 1) 

X,Y sampling intervals, m (see fig. 1) 

y equivalent ins tantaneous field of view (IF OV) , m (see fig 

6 (x) unit impulse function, m" 1 

S(x,y) two-dimensional unit impulse function, m" 2 
0 angle in the xy-plane, deg 

y (expected) mean value, Wm" 2 

0,0) spatial frequency along x- and y-axis, respectively, m"* 1 

p = \\j2 + 0)2 j 1 /2 

a standard deviation, Wm" 2 

a standard deviation of smoothed measurements 

T(x,y) spatial response (or point-spread function) 


. 3) 


(cycles/m) 
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< H O <0 


(U,oj) spatial frequency response (or modulation transfer function (MTF) ) 

(•) autooovariance, W^nT^ 

(•) Wiener (or power density) spectrum, W 2 m -2 (w 2 m - Vcycles 2 m -2 ) 

II(*) rectangle function 

111 (•) sampling or comb function 

* convolution 

Subscripts: 

a aliased noise, or sideband component generated by sampling 

e electronic filter 

i. lens 

L radiance fluctuations 

m measurement 

p photosensor aperture 

r interval between radiance fluctuations 

R continuous signal reconstructed from discrete measurements 

s sufficiently sampled signal components 

DEFINITIONS, ASSUMPTIONS, AND FORMULATIONS 
Measurement Process 

This paper analyzes the linear, space- invariant measurement process shown 
in figure 1 that transfers the (continuous) radiance fluctuations L(x,y) into 
a (discrete) signal m(x,y;X,Y) as defined by the equations 

m (x,y) = L(x,y) * T(x,y) (la) 


m (x,y; X,Y) 


m (x,y) 



CL(x,y) 


* T(x,y) ] 



(lb) 
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where the symbol * denotes convolution; T(x,y) r the spatial response (or 
point-spread function) of the radiometer; and 1 1 1 ' the sam P lin 9 (ref. 11) 

or comb (ref. 12) function given by 


JLLL(r»“) = IEZ 4z ■ *' “ " n ) = ™ IEZ SI «(x - xm, y - Yn) 


VX Y i 


m=-°° n=-°° 


m=-°° ns-" 


The discrete measurements m(x,y;X,Y) are an infinite array of unit impulse 
functions on the rectangular lattice with spacing X and Y in the x- and 
y-directions, respectively. Rectangular (Cartesian) coordinates (x r y) in the 
object or photosensor aperture plane are used as the reference for the measure- 
ment coordinates. 

Salient properties of the measurements m(x,y;X,Y) defined by equa- 
tion (1b) are often more conveniently evaluated in the frequency domain than 
in the spatial domain. The functions g(U,u)) and g(x,y) are denoted as the 
Fourier transform pairs given by 


w 

g(U,w) = g (x,y)e“-*-2Tr(Ux+0)y) ^ 


(2a) 


g(x,y) = g(V,w)e*2 7T ( x ' J+ y' t °) <JU do) 


(2b) 


The Fourier transform of equation (1b) yields 


(U,u); X, Y) = L(U,OJ) T(U,to) * XY j 1 1 (XU,Yco) 


m 


(3a) 


where 


oo oo 


^ ^ _2_ _!L_ / m n^ 

XY 1 1 1 (XU,Y 0 J) = XY > 2 <$(XU - Yo) - n) = 2 > P , 10 


m»-" n=-" 


m=-" n=-" 


X Y J 
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Equation (3a) can be equivalently written as 


m(V),u); X,Y) 


00 oo 





m 

0) 

X 




(3b) 


Radiance Fluctuation 

The magnitude of aliased noise in continuous reconstructions of discrete 
radiometric measurements depends on the properties of the radiance fluctuations 
L(x,y) as well as on the spatial response T(x,y) and the sampling intervals 
(X,Y) . In this paper a random radiance fluctuation that is both homogeneous 
and isotropic (ref. 13) is assumed, so that the variance of L(x,y) is inde- 
pendent of (x,y) , and the autocovariance (and autocorrelation) is a function 
only of the relative distance 


r = [(x-, - x 2 ) 2 + (yi - y 2 J 2 ] 1 / 2 


(4) 


between points (x-| ,y-| ) and (x 2 ,y 2 ) . It follows that the assumed distribution 
of L(x,y) is circularly symmetric. 

Furthermore, the autocovariance of L(x,y) is assumed to be 


* L ( r > 


Or . 2 


,-r/li r 


(5) 


Then, the associated Wiener (or power density) spectrum, which is the circu- 
larly symmetric Fourier transform (i.e. , the Hankel transform) of $£,(r), is 


$ L (u,co) = $ L (p) 


2iry r 2 a L 2 


[1 + (2trvi r p) 2] 3/2 


( 6 ) 


where p 2 = U 2 + w 2 . (See fig. 2.) Equations (5) and (6) can be derived by 
assuming that L(x,y) is a set of pulses whose random location in the xy-plane 
obeys the Poisson probability density function with the (expected) mean width 
U r and whose random magnitude obeys the Gaussian probability density function 
with the (expected) mean and variance and a L *, respectively. (For deri- 

vations and further details, see refs. 14 and 15; for a detailed one-dimensional 
analysis, see ref. 16.) The variance a L 2 and the Wiener spectrum $ L (U,a)) 
are related by 
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(7) 


a L 2 = $ L (0) 


00 



—oo 


$ L (U,0)) 


du dw = 2TI 



p$ L (p) dp 


That is, dj, 2 is the value of the integrated Wiener spectrum. 

Radiometer Response 

The spatial frequency response (or modulation transfer function (MTF) ) 
t(u,o)) of a radiometer can generally be characterized as 


T(U,U) = T£(U,10) Tp(U,lO) T e ((jj) (8) 


A A a 

where T£(u,o)), Tp(U,o>) , and T e (a)), respectively, are the frequency response 
of the objective lens, photosensor aperture, and electronic filter; ti) is the 
spatial frequency component in the line-scan direction. The photosensor aper- 
ture and its distance from the lens define the IFOV of the radiometer. The 
lens aperture and IFOV govern the total amount of radiance focused onto the 
photosensor aperture; however, unlike for high-resolution imaging devices, the 
(diffraction-limited) spatial frequency response of the lens usually does not 
significantly affect the frequency response of relatively low-resolution radi- 
ometers and is neglected in this paper. 

Figure 3 illustrates salient characteristics of three photosensor aperture 
shapes. These apertures have equivalent areas, namely Try^/4. Thus for each 
aperture. 


and 



Tp(x,y) dx dy = T p (0,0) 


1 



((x,y) aperture area) 


(elsewhere) 


It is convenient to regard the diameter y of the circular aperture as the 
equivalent IFOV of the three apertures. The spatial frequency response 
Tp(U,oo) of the aperture shapes illustrated in figure 3 is given by the 
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following expressions , in which the shape parameters a and b are the x- and 
y-axis intercepts, respectively, of the apertures, and sine x = sin (ttx)/itx: 

(a) Circular - 


J-j (2TTap) 

t p (U,cd) = T p (p) = 

F v Trap 


(9) 


where a 


Y 

2 


(b) Cosine - 


T p (U,U)) 


4 pTT/2 

— I cos (aU0) sine (2bto cos 2 9) cos 2 9 d9 
tt J 0 


( 10 ) 


TT 

where a = — y 

4 


b 


Y 

2 


(c) Diamond - 


Tp(U,u)) = sine (aU + boo) sine (aU - b 03 ) 


( 11 ) 


TT 

where a = — Y 

4 


b 


Y 

2 


The effect of various photosensor aperture shapes and profiles on the 
image quality obtained with line-scan devices has been investigated most nota- 
bly by Mertz and Gray (ref. 1) for facsimile scanners and by Schade (ref. 2) 
for television cameras. One of their conclusions was that all reasonable pro- 
files and shapes of equivalent size result in about equal blurring, but that 
some of the profiles and shapes tend to suppress aliasing better than others. 
Katzberg et al. (ref. 6) demonstrated that the cosine and diamond shapes are 
both superior to the circular shape in suppressing aliasing. 

The cosine and diamond apertures are shaped to suppress side lobes along 
the U-spatial frequency axis normal to the line-scan direction. The exact 
aperture shape along the line-scan direction is less important in continuously 
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scanning devices, since the effective w-spatial frequency response can be 
reshaped by an electronic low-pass filter. The function given by 


T e <(0) 



(M < i/y) 

( 12 ) 

(|w| > I/y) 


and shown in figure 4 was chosen to represent the spatial frequency response 
of an electronic filter. The cutoff frequency 1 /y provides a useful spatial 
frequency passband for the radiometer that is approximately circularly symmet- 
ric. In practice, it might often be preferable to eliminate aliasing along the 
line-scan direction by selecting a cutoff frequency equal to 1/2Y. 


Signal and Aliased Noise 

The discrete measurements m(x,y;X,Y) defined by equation (lb) have 

/% 

a mean value y L and a Wiener spectrum 4> m (u,a);X,Y) given by 


$ m (U,w;X,Y) = ® S (U,U» + $ a (U,U);X,Y) 


(13) 


The term $ s (u,w) represents the (desired) fundamental signal component cen- 
tered at the location (U,u)) = (0,0) and is given by 


$s(U,oj) = $ L (U,u» |x(u,w) | 2 


( 14 ) 


The term $ a (u,oj;X,Y) represents the (undesired) sampling sideband components 

/ m n\ 

centered at the locations (U,w) = — , (m,n) / (0,0), and is given by 

\x y/ 


” / m n> 

$ a (U,U;X,Y) = > > ® s O - -, W - - 

m=-°° n=-°° \ * ^ / 


(15) 


(m,n) / (0,0) 


The development of equation (15) from equations (1) is standard. (See for 
example, refs. 10 and 13.) 
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The frequency passband of the sampling lattice with spacing (X,Y) , or 
briefly the sampling passband, is given by the rectangle function (ref. 11) 


n (xo,yoj) = 


r 1 1 \ 

M < — , | to | < — 

, 2X 2Y/ 


( 16 ) 


(elsewhere) 


a g 2 of the fundamental signal component within the sampling 


^ pl/2X pl/2Y ^ 

$ s (U,w) n(XU,Y(D) du do) = \ \ $ S (U, 0 )) dv) dw (17) 

J-1/2X -1/2Y 

and the variance a a ^ of the sampling sideband components within the sampling 
passband is 

r l/2X r l/2Y A 

$ a (U,0);X,Y) n(XU,YCO) dU d (0 = \ \ $ a (U,U);X,Y) dl) dU) (18) 

J -1/2X J-1/2Y 


If the Wiener spectrum $ L (u,to) of the radiance fluctuations or the 

A 

spatial frequency response t(U,(jl)) of the radiometer is contained within the 

A 

sampling passband (i.e., if $ s (t) f aj) = 0 outside the sampling passband), 
then a a 2 = 0, and the signal is said to be sufficiently sampled. Otherwise, 
the signal is insufficiently sampled, and displaced frequency components from 
the sidebands fall into the sampling passband, as illustrated in figure 5. 

These displaced frequency components cannot be distinguished in practice from 
the desired signal - they tend to behave, and are treated here, as additive 
noise. This treatment follows Blackman and Tukey (ref. 17), who define aliased 
noise as those components of the signal that are contained in the signal only 
when insufficient sampling occurs. Insufficient sampling occurs frequently in 
practice because the spatial frequency response t(u,u)) of optical systems 
tends to decrease very gradually and to oscillate around zero. 

The function m R (x,y) is denoted as a continuous reconstruction of the 
discrete measurements m(x,y;X,Y) ; the spatial resolution of this reconstruc- 
tion is defined as the resolution cell (resel) R. If the rectangle function 
II(XU,YW) given by equation (16) is the reconstruction filter, then m R (x ,Y) 
is the inverse Fourier transform of m(u,to;X,Y) II(XU,Yo); that is. 



The variance 
passband is 
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oo oo 

»R(x,y) = m(mX,nY) sine 

ms-oo n =-°° 


where ra(mX,nY) is the value of m(x,y) given by equation (la) at the sam- 
pling points (mX,nY). The area of the corresponding resel R is XY. The 
Shannon-Whittaker sampling theorem states that the reconstruction m R (x,y) is 
equal to the signal m(x,y) if m(x,y) has been sufficiently sampled (i.e. , 
a a 2 = 0) . 

It is occasionally desirable to produce spatially smoothed reconstructions, 
for example, to gain increased radiometric accuracy of the reconstructed signal 
at the cost of spatial resolution. Spatial smoothing can be accomplished by 
narrowing the passband of the reconstruction filter. It is convenient to 
express the reconstruction filter as a function of the equivalent IFOV y as 



II(yMU,yNa)) = 


[M <— , M < — 

l 2KM 2KH, 


(20) 


(elsewhere) 


where M and N are constants. The smoothed reconstruction is 


mR(x,y) 


XY 

yMyN 


CO oo 

> > m(mX,nY) 

m=-°° n=-°° 


sine 




(21) 


The area of the corresponding resel R is y 2 MN. The variance of the signal 
and aliased noise of %(x,y) is given by 


r l/2yM r l/2yN A 

o s 2 - \ \ $ s (u,u>) du did 

J -l/2yM ^-l/2yN 


and 


r 1/2yM pl/2yN „ 

a a 2 = \ \ $ a (u,w ; x,Y) du du 

J -l/2yM ^-1/2yN 


( 22 ) 


(23) 
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COMPUTATIONAL RESULTS 


Figures 6 to 8 characterize the (normalized) standard deviations of the 
signal, aliased noise, and signal- to-noise ratio of spatially unsmoothed and 
smoothed reconstructions of the discrete measurements m(x,y;X,Y) . The stan- 
dard deviations of the signal 0 S and the aliased noise a a are normalized 
with respect to the standard deviation cr^ of the radiance fluctuations; the 
sampling intervals X and Y and the mean spatial width y r of the radiance 
fluctuations are normalized with respect to the equivalent IFOV y. These 
normalizations permit the computational results to be applied to a wide range 
of practical situations. 

Figure 6 presents the variation of the signal, aliased noise, and signal- 
to-noise ratio versus sampling intervals for various mean spatial widths of the 
radiance fluctuations. These results characterize the unsmoothed signal recon- 
structions depicted by equations (16) to (19). 

Figure 7 presents these variations for the smoothed signal reconstructions 
formulated by equations (20) to (23) . It may be noted that the conditions (and 
results) for X/y = Y/y « 1 and M = N = 1 are the same as those in figure 6 
for X/y = Y/y = 1 . 

Figure 8 presents the same computational results as figure 7. The differ- 
ence is that figure 7 shows variations in performance plotted against sampling 
interval while figure 8 shows variations plotted against mean spatial width of 
the radiance fluctuations. 


CONCLUSIONS 

1 . The accuracy with which sufficiently sampled components of random radi- 
ance fluctuations (or fields) can be reconstructed, as accounted for by 0 s /0l, 
tends to be independent of reasonable photosensor aperture shapes of equal area; 
however, the degradation of this accuracy by insufficiently sampled components 
of the radiance fluctuations, as accounted for by the aliased noise tfa/°L' 
depends strongly on these aperture shapes. This conclusion is consistent with 
the observation by Mertz and Gray (Bell Syst. Tech. J., vol. 13, July 1934) and 
by Schade (J. SMPTE, vol. 64, Nov. 1955) that all reasonable photosensor aper- 
ture shapes and profiles of equivalent size result in about equal blurring, but 
that some shapes and profiles tend to suppress aliasing better than others. 

2. Spatial smoothing (i.e., low-pass filtering) of discrete measurements 
decreases both the accuracy with which the sufficiently sampled radiance fluc- 
tuations (i.e., ^g/Ck) can be reconstructed in the absence of aliased noise 
and the degradation of this accuracy by the presence of aliased noise (i.e., 
a a /a L ) . There is a net gain in the signal- to-noise ratio ag/a a of 
smoothed reconstructions because a a tends to decrease more rapidly with 
smoothing than a s . 

3. Aliased noise is a large source of error in continuous reconstructions 
of discrete measurements obtained with contiguous coverage, unless the mean 
spatial width of the (random) radiance fluctuations is more than an order of 
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magnitude larger than the instantaneous field of view (IFOV) of the radiometer 
(i.e., unless n r /y >10). For example, for a circular IFOV and for unsmoothed 
reconstructions, the signal- to- noise ratio o s /a a is only 16 for M r /y =9, 

5 for y r /y = 1, and 1.4 for y r /y m 1/9. 

4. A decrease of about 2 orders of magnitude in aliased noise and a corre- 
sponding increase in signal-to-noise ratio are attained by modestly decreasing 
the sampling interval (by a factor of 0.7) and by shaping the effective IFOV 
(i.e., photosensor aperture and electronic filter) of a scanning radiometer. 
This improvement tends to be independent of the statistical properties of the 
radiance fluctuations. For example, the signal-to-noise ratio o s /o a for 
unsmoothed reconstructions is increased from 16 to 514 for li r /y = 9, from 5 

to 139 for 1J r /y = 1, and from 1.4 to 26 for y r /y = 1/9. 

5. A much more modest improvement is attained at the cost of spatial reso- 
lution by smoothing (i.e., low^pass filtering) of the discrete measurements . 

The signal-to-noise ratio o s /a a obtained with a circular IFOV tends to 
increase by about the same factor by which the reconstructed resolution is 
decreased. The diamond IFOV (with electronic filter) provides a more favorable 
trade-off. 

6. For unsmoothed reconstructions, the aliased noise Oa/ a L decreases and 

the signal-to-noise ratio cr s /a a increases with decreasing sampling interval. 
The rate of change varies: relatively rapid improvements can be attained by 

decreasing the sampling intervals X/y = Y/y from 1 to 0.7; further decreases 
in sampling intervals provide less improvement. For smoothed reconstructions, 
the variation of ^a/°h 3°^ o s /o a with sampling intervals can exhibit 

intermediate dips and peaks, respectively. This occurs when the zero crossing 
of the neighboring sampling sidebands falls into the (narrow) low-pass recon- 
struction filter. 

7. The standard deviations of the signal Og/a ^ and the signal-to-noise 
ratio o s /cj a decrease monotonically as the mean spatial detail y r decreases, 
whereas the standard deviation of the aliased noise Og/o-^ reaches a broad 
maximum when the mean spatial detail \i x is slightly smaller than the IFOV y 
(i.e., when 0.1 < y r /y < 1). 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 13, 1980 
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(a) CcmiDon measurement process: contiguous coverage (i.e., x/y = y/y = 1) 

with discretely stepped, circular IFOV. 

Figure 5.- Measurement characteristics in spatial frequency domain. Only 
neighboring sampling sidebands shown here have been included in compu- 
tation of aliased noise; more accurate computations that include addi- 
tional sidebands would yield slightly larger magnitudes of aliased noise. 
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(a) Circular IFOV without *l*ctronic filtor. (b) Cosin* IFOV with *l*ctronic filter. (c) Diamond IFOV with *l*ctronic filter. 


Figure 6.- Variation of standard deviation of normalized signal tfg/a L , aliased noise tf a /a L , and 

signal-to-noise ratio cr s /a a , versus sampling intervals X/y - Y/y for various mean spatial detail 
widths y r /y of radiance fluctuations L(x,y)* These results characterize unsmoothed reconstruc- 
tions formulated by equations (16) to (19). Area of resel R is XY. 





X/ r =Y/ y X/ X = t/ X X/ X = Y/y X/ X =Y/ r 

|b) Diamond IFOV (with lUcIronic lill.r) 

Figure 7.- Variation of standard deviation of normalized aliased noise cr a /a L , and of signal-to-noise 
ratio 5 g /cr a versus sampling intervals X/y = Y/y for various mean spatial detail widths 
y r /y of radiance fluctuations L(x,y) . These results characterize smoothed reconstruction formu- 
lated by equations (20) to (23) . Area of resel R is y 2 MN; number of measurements encompassed 
by R is y 2 MN/XY (e.g., MN if X/y = Y/y =1, and 4MN if X/y = Y/y = 0.5). 











Figure 8.- Variation of standard_deviation of normalized signal a s /a L and aliased noise 0 a /a L , and 
of signal-to-noise ratio a s /5 a versus mean spatial detail widths M r /Y* These results char- 
acterize smoothed reconstructions formulated by equations (20) to (23) . 
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